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mators  of  the  failure  rate  function  and  survival  probability  were  developed 
under  the  assumption  of  Increasing  failure  rate  using  both  maximum  likelihood 
and  Bayesian  approaches.  These  particular  results  have  attracted  wide  atten¬ 
tion  due  to  their  generality  and  applicability  in  survival  analysis  and 
reliability  estimation  from  arbitrarily  right-censored  data.  Also,  consis¬ 
tency  results  for  both  univariate  and  multivariate  kernel  estimates  for 
probability  density  functions  and  regression  functions  were  obtained  using 
techniques  and  results  of  function-space  probability  theory.  Stochastic 
convergence  results  for  weighted  sums  of  random  elements  in  various  function 
spaces  were  also  obtained.  Sequential  procedures  were  developed  and  analyzed 
which  provided  interval  estimators  of  the  parameter  of  interest  after  testing 
certain  hypotheses.  Various  robustness  and  nonparametric  methods  for  incom¬ 
plete  samples  and  broken  samples  were  studied,  as  well  as  nonparametric 
repeated  significance  testing  and  randomized  block  designs.  Thus,  maintenanoej 
policies  and  development  of  new,  more  reliable ,  equipment  may  be  formulated  | 
using  the  statistical  procedures  and  theory  from  these  results. 
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I.  Introduction 

This  document  reports  in  detail  the  work  performed  and  other  research 
activities  of  the  principal  investigators  during  the  period  from  June  1,  1979, 
to  May  31,  1981,  under  contract  number  ^■BR-F49620-79-C-0140. 

In  Section  2  a  comprehensive  statement  of  the  research  objectives 
during  this  period  is  given,  and  in  Section  3  a  detailed  report  of  the  accom¬ 
plishments  and  highlights  of  the  significant  results  obtained  is  presented. 
Section  4  contains  a  complete  listing  of  approximately  thirty-one  publications 
which  have  resulted  from  these  research  efforts.  The  personnel  supported  under 
this  contract  are  listed  in  Section  3,  and  a  list  of  the  professional  interac¬ 
tions — talks  at  conferences  and  meetings,  colloquia  and  seminars,  etc. —  is 
presented  in  Section  6.  Section  7  gives  an  indication  of  the  significance  of 
the  results  to  military  application.  Finally,  other  professional  activities 
of  the  principal  investigators  are  outlined  in  Section  8. 


Approved  for  publlo  release  f 
distribution  unlimited. 


2 


Research  Objectives  During  the  Contract  Period 

The  research  objectives  of  this  contract  can  be  divided  into  five 
broad  catagories  or  areas: 

(A)  Survival  analysis  and  life  testing; 

(B)  Probability  density  estimates  and  laws  of  large  numbers; 

(C)  Estimation  after  testing; 

(D)  Robustness  and  distribution-free  procedures; 

(E)  Stochastic  systems  theory. 

The  research  problems  which  were  considered  under  this  contract  in  these  five 
areas  will  be  outlined  in  this  section.  The  progress  toward  their  solution 
and  significant  results  will  be  described  in  Section  3  of  this  report. 

A.  Survival  Analysis  and  Life  Testing 


The  main  thrust  of  the  research  in  this  area  was  to  develop  estimators 
of  survival  probability  or  reliability  in  nonparametric  settings  or  in  a  non- 
parametric  Bayesian  framework  which  w.  re  much  better  than  existing  estimators 
in  the  sense  of  their  smoothness  and  convergence  properties.  The  main  objective 
in  this  area  was  to  develop  such  estimators  for  arbitrarily  right-censored 
observations  under  no  distributional  assumptions  other  than  assumptions  con¬ 
cerning  the  failure  rate.  In  particular,  since  increasing  failure  rate  (IFR) 
distributions  are  very  commonly  used  to  model  wear-out  failures,  the  distribu¬ 
tional  assumption  was  simply  that  the  underlying  life  time  distribution  was 
IFR.  One  of  the  most  promising  approaches  to  these  problems  was  to  consider 
(nonparametric)  maximum  likelihood  estimation  of  an  increasing  failure  rate 
function  r(t),  t  2  0,  based  on  arbitrarily  right-censored  observations  from 


the  life  distribution  F  and  then  to 
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F(t)  *  1  -  F(t)  -  exp(-  /q  r(u)du)  from  the  maximum  likelihood  estimate  of  r. 
The  other  promising  avenue  of  attack  for  these  problems  was  to  assume  a  prior 
distribution  on  the  space  of  all  increasing  failure  rate  functions  and  obtain 
a  nonparametrlc  Bayes  estimator  of  F  or  F  which  was  a  continuous  function.  A 
type  of  Poisson  "shock  model"  was  successfully  used  as  a  prior  distribution. 

Another  objective  in  this  area  was  to  develop  powerful  nonparametric 
tests  for  symmetry  of  a  bivariate  distribution  function  with  observations 
which  are  subject  to  arbitrary  right  censorship.  The  usual  sign  test  is  not 
the  best  test  to  use  in  the  case  of  censoring  since  a  pair  of  observations,  one 
of  which  is  censored,  may  not  be  comparable. 

Also,  under  the  parametric  inverse  Gaussian  and  Birnbaum-Saunders 
models,  it  was  desired  to  obtain  Bayes  estimators  of  the  reliability  function. 
Approximate  prediction  intervals  for  the  mean  of  future  lifetimes  from  the 
inverse  Gaussian  model  were  to  be  investigated,  as  well  as  the  problem  of 
estimating  reliability  under  growth  assumptions  for  the  two-parameter  gamma 
distribut ion. 

B.  Probability  Density  Estimates  and  Laws  of  Large  Numbers 

This  area  of  the  research  concentrated  on  obtaining  consistency  of 
probability  density  estimates  and  developing  laws  of  large  numbers  which 
would  have  applications  in  the  estimation  of  probability  density  functions. 

In  particular,  the  kernel  density  estimation  techniques  were  examined  with 
the  research  objective  of  obtaining  a  wider  class  of  kernel  functions  and 
more  general  bandwidth  sequences  which  could  be  used  while  preserving  the 
stochastic  convergence  of  the  estimate  to  the  true  underlying  probability 
density  function.  The  development  of  multivariate  kernel  density  estimates 
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with  convergence  properties  which  were  similar  to  the  univariate  estimates 
was  also  a  research  objective. 

In  the  kernel  density  estimation,  a  kernel  function  (usually  a 
probability  density  function  which  the  experimenter  may  choose)  is  centered 
at  each  observed  sample  point  and  is  scaled,  depending  on  the  sample  size,  so 
that  its  spread  decreases  and  becomes  more  concentrated  around  the  sample 
observation  as  the  sample  size  increases.  The  estimate  is  then  an  average  of 
these  kernel  functions  which  are  centered  at  the  sample  observations.  Unfor¬ 
tunately,  the  changing  bandwidth  (as  the  sample  size  increases)  forces  consi¬ 
deration  of  a  Marcinkiewicz-Zygmund  type  of  such  an  average.  Hence,  laws  of 
large  numbers  in  function  spaces  containing  the  possible  kernel  functions 
would  yield  consistency  of  these  estimates  and  was  a  major  research  objective 
of  the  grant.  Since  the  scaled  kernel  functions  changed  as  the  sample  size 
changed,  the  appropriate  formulation  involved  weighted  sums  of  arrays  of 
funct ion-val ued  random  variables.  Also,  the  possible  different  measures  of 
dispersion  of  the  kernel  estimate  f^(t)  from  the  probability  density  function 
f(t)  necessitated  the  consideration  of  laws  of  large  numbers  for  Banach  spaces 
and  the  function  space  D[0,1]  as  a  research  objective. 

C.  Estimation  After  Testing 

Suppose  F  is  a  distribution  function  which  is  stochastically  nonde- 
creasing  in  6.  It  was  desired  to  obtain  various  sequential  procedures  which 

allow  quick  rejection  of  the  null  hypothesis  H  :  6  >  0  when  the  alternative 

'  o  o 

H^:  0 f  0 1  (0^  >  Gj)  is  true.  In  addition,  these  procedures  would  provide  an 
accurate  interval  estimate  of  the  parameter  0  when  the  null  hypothesis  is 
accepted.  The  results  would  hold  in  particular  for  an  exponential  distri¬ 


bution  with  mean  0. 
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D.  Robustness  and  Distribution-free  Procedures 

Several  of  the  research  objectives  were  in  this  category.  One  is 
concerned  with  the  distributional  robustness  when  there  are  missing  obser¬ 
vations  in  the  samples.  The  second  one  was  to  establish  the  robustness 
(asymptotic  conservativeness)  of  a  K-sample  Kruskal-Wallis  test  under  certain 
departures  from  mutual  independence  of  K  univariate  samples.  Another  objective 
was  to  Investigate  an  asymptotically  distribution-free  simultaneous  confidence 
region  of  treatment  differences  in  a  randomized  complete  block  design.  Also, 
it  was  desired  to  obtain  nonparametric  repeated  significance  tests  for  two- 
sample  survival  data  to  test  equality  of  the  two  distributions  and  to  estimate 
location  difference  for  fragmentary  samples.  The  estimation  of  the  ratio  of 
the  scale  parameters  of  two  distributions  based  on  arbitrarily  right -censored 
observations  from  each  distribution  was  also  to  be  considered  in  detail.  Very 
good  results  were  obtained  in  all  of  these  areas  of  investigation. 

E.  Stochastic  Systems  Theory 

The  main  objective  in  the  study  of  stochastic  systems  was  to  obtain  very 
general  conditions  for  the  existence  and  stability  of  solutions  to  a  class  of 
stochastic  systems  which  can  be  described  by  stochastic  integral  equations 
involving  McShane  integrals.  This  class  includes  stochastic  integro-dif ferent ial 
systems.  Frequency  domain  conditions  for  the  stability  of  the  solution  were  to 


be  considered. 
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3.  Results  of  the  Research  Effort 

In  this  section,  a  substantive  summary  of  the  significant  accomplishments 
towards  achieving  the  research  objectives  of  this  contract  as  outlined  in  Section  2 
of  this  report  will  be  given.  The  specific  research  papers  containing  the  results 
will  be  listed  in  Section  4.  Copies  of  these  research  papers  have  been  forwarded 
to  the  Program  Manager  as  they  were  submitted  for  publication. 

A.  Survival  Analysis  and  Life  Testing  / 

One  of  the  most  significant  results  in  this  area  is  paper  number  1  in  the 
listing  of  Section  4.  This  paper  had  received  wide  attention  from  workers  in  this 
area,  even  before  it  had  appeared  in  print,  since  it  was  presented  at  the  annual 
Institute  of  Mathematical  Statistics  Meeting  in  Washinton,  D.C.,  in  August,  1979. 
The  result  concerns  the  (nonparametric)  maximum  likelihood  estimation  of  an  in¬ 
creasing  failure  rate  function  and  the  corresponding  distribution  function  based 

; 

on  a  set  of  observations  subject  to  arbitrary  right  censorship.  The  estimator  is 
developed  using  the  techniques  of  isotonic  regression,  is  defined  everywhere  on 
the  positive  real  line,  is  a  smooth  function,  and  decreases  to  zero  as  time  increaj 
The  small  sample  properties  of  the  estimator  were  compared  with  the  Kaplan-Meier 
estimator  by  a  Monte  Carlo  study  for  Weibull  distributions,  and  the  results  indi¬ 
cated  that  the  maximum  likelihood  estimator  is  superior  to  the  Kaplan-Meier  esti¬ 
mator.  This  estimator  is  applicable  to  a  broad  spectrum  of  problems  due  to  the 
generality  of  the  censoring  mechanism.  For  example,  medical  trial  data  as  well  as 
industrial  life  test  data  may  be  analysed  with  this  estimation  procedure.  In 
particular,  since  no  assumptions  are  required  concerning  the  exact  lifetime  dis¬ 
tribution,  this  estimator  is  especially  applicable  to  the  assessment  of  military 
equipment  reliability  and  maintenance.  Such  equipment  which  is  subject  to  wear- 
out  may  be  assumed  to  follow  an  1FR  life  distribution,  and  this  is  the  only 
assumption  necessary  for  the  estimation  procedure  developed  here. 


A  second  significant  result  in  this  area  is  contained  in  paper  number  5 
in  Section  4.  Bayesian  nonparametric  estimators  of  the  survival  function,  the 
failure  rate  function,  and  the  density  function  are  obtained  using  Poisson 
jump  processes  to  generate  prior  probability  measures  on  the  space  of  increasing 
failure  rate  functions.  The  jump  processes  are  intuitively  appealing  and  have 
a  meaningful  physical  interpretation  as  a  shock  model,  where  "shocks"  occurring 
to  the  system  increase  the  failure  rate  by  a  positive  amount  at  the  random 
instant  of  occurrence.  The  Bayesian  nonparametric  estimators  are  also  easily 
derived  for  these  priors  for  arbitrarily  right-censored  observations.  These 
estimators  are  continuous  and  are  thus  more  appealing  in  this  sense  than 
estimators  of  other  authors  who  used  Dirichlet  processes  as  prior  distributions 
for  (discrete  with  probability  one)  distributions.  The  same  comments  made  on 
the  applicability  of  the  estimator  in  the  last  paragraph  can  be  made  for  these 
nonparametric  Bayes  estimators.  Also,  along  these  same  lines  when  the  Dirichlet 
process  is  not  completely  specified,  several  sequences  of  nonparametric  empir¬ 
ical  Bayes  estimators  of  survival  probability  were  considered  (paper  number 
18  in  Section  4)  and  their  asymptotic  optimality  with  respect  to  a  Dirichlet 
process  was  investigated. 

In  testing  a  bivariate  life  distribution  for  symmetry  based  on  arbitrar¬ 
ily  right-censored  observations,  a  conditionally  distribution-free  test  was 
developed (paper  number  20  in  Section  4).  A  permutation  test  using  Gehan's 
scores  is  proposed  which  is  much  more  powerful  than  the  sign  test.  The  power 
of  the  test  was  compared  with  that  of  the  sign  test  by  computer  simulations 
using  the  Marshal l-01kin  bivariate  exponential  model. 

For  the  inverse  Gaussian  distribution  as  a  parametric  lifetime  model, 

Bayes  estimators  of  the  reliability  function  have  been  obtained  (paper  number 
12  in  Section  4).  For  the  case  that  the  mean  lifetime  is  known,  Bayes  estimators 
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of  the  reliability  are  obtained  with  respect  to  squared  error  loss  for  Jeffreys' 
non in formative  prior  and  for  the  gamma  family  (natural  conjugate  family)  of 
priors  for  the  scale  parameter.  In  the  case  that  both  the  mean  and  scale  para¬ 
meter  are  unknown,  an  estimator  of  reliability  is  suggested  which  is  based  on 
the  Bayes  estimator  for  the  first  case.  This  estimator  compares  favorably 
with  the  maximum  likelihood  and  minimum  variance  unbiased  estimators. 

The  Birnbaum-Saunders  fatigue  life  distribution  with  shape  parameter 
a  and  scale  parameter  g  was  also  considered  (paper  number  19  in  Section  4). 

The  scale  parameter  g  is  also  the  median  lifetime,  and  assuming  that  g  was 
known,  Bayes  estimators  of  the  reliability  function  were  obtained  for  a 

family  of  proper  conjugate  priors  as  well  as  for  Jeffreys'  vague  prior  for  a- 

* 

When  both  a  and  g  were  unknown,  a  modified  Bayes  estimator  Rg(t)  of  the 
reliability  function  was  proposed  using  a  moment  estimator  of  g.  In  addition 
to  being  computationally  simpler  than  the  maximum  likelihood  estimator  of 

the  reliability  function,  Monte  Carlo  simulations  for  small  samples  showed 

* 

that  Rg(t)  was  better  than  the  method  of  moments  estimator  for  all  a  and  as 
good  as  the  maximum  likelihood  estimator  for  small  a  in  the  sense  of  mean 
square  errors. 

The  problem  of  predicting,  on  the  basis  of  an  observed  sample  from 
an  inverse  Gaussian  distribution,  the  mean  of  a  future  random  sample  from 
the  same  distribution  was  also  considered  (paper  number  13  in  Section  4). 
Approximate  prediction  intervals  were  developed  since  it  seemed  to  be  impos¬ 
sible  to  derive  exact  intervals.  The  accuracy  of  the  approximate  intervals 
was  investigated  by  extensive  Monte  Carlo  simulations  and  found  to  be  very 
good.  These  results  are  useful  for  predicting  the  next  first  passage  time 
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for  a  Brownian  motion  with  positive  drift  or  the  failure  time  of  an  item 
having  inverse  Gaussian  life  distribution. 

Finally,  in  this  area  the  estimation  of  the  reliability  function 
under  development  testing  was  investigated  (paper  number  22  in  Section  4). 

Changes  are  often  made  at  various  stages  during  the  development  of  a  system 
to  correct  design  weaknesses.  This  results  in  reliability  growth,  an  increase 
in  reliability  and  mean  life  of  such  systems  as  development  continues.  Assuming 
that  the  failure  distribution  at  stage  i  of  K  development  stages  is  gamma  with 
parameters  and  g^,  approximate  maximum  likelihood  estimates  of  c^,  g^,  i*l,...JC, 
were  obtained  subject  to  the  conditions  and  g  <  . . .  <g  ,  which  imply 

reliability  growth.  An  iterative  procedure  involving  two  constrained  nonlinear 
optimization  problems  was  presented  for  obtaining  the  approximate  maximum  like¬ 
lihood  estimates.  The  solutions  to  the  constrained  optimization  problems  were 
relatively  simple  to  compute  using  techniques  of  isotonic  regression,  and  the 
iterations  tended  to  converge  rather  quickly. 


Probability  Density  Estimators  and  Laws  of  Large  Numbers 


In  paper  number  2  in  the  listing  of  Section  4,  the  uniform  consistency 


of  the  multivariate  kernel  density  estimate,  f^(Xj , . . . ,x^) ,  was  obtained  by 
showing  that  | | f^(x^ , . . . ,x^)  -  f (x ^ , . . . ,x^) j | m  converges  completely  to  zero  for 
a  large  class  of  weight  functions  and  under  mild  conditions  on  the  bandwidth 


sequences  b^(n)'s.  This  significant  result  provided  for  the  uniform  consistency 
~  hn(xj » •  •  •  > 

of  the  estimate  m  (x  ,...,x.)  =  — — — - for  the  regression  function 

"  . V 

★ 

m(x1,...,xk)  *  E(Y | X j  =  xx . ^  =  x^)  where  f(xj . x^)  =  / f  (xj , . . . ,xk,y) dv 

and  hCXj,...,^)  =  Jyf  (x j , . .  . , x^,y) dy.  In  these  results  the  rates  of  conver¬ 


gence  need  not  be  uniform  in  each  coordinate  and  the  broad  class  of  possible 
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kernels  included  Epanec hnikov ' s  optimal  kernel  and  was  characterized  by  a 
smoothness  condition. 

In  papers  number  24  and  25  in  the  listing  of  Section  4,  the  stochastic 
convergence  of  weighted  sums  of  arrays  of  random  elements  in  Banach  spaces  was 
obtained  under  varying  conditions.  One  striking  result  was  that  the  convergence 
of  these  weighted  sums  was  necessary  and  sufficient  for  the  geometric  property, 
type  P»  °f  the  Banach  space.  As  corollaries,  these  results  have  Banach  space 
versions  of  several  random  variables  results.  In  paper  number  9  it  is  shown 
in  the  application  of  these  results  that  the  consistency  in  the  L^-norm  of 

ZOO  2  —  q 

(b  n)  <  00  for  some 

i  n 

n=  1 

q  >  1  where  n  is  the  sample  size  and  b^  is  the  bandwidth.  The  L^-consist ency 

for  a  broad  class  of  kernel  estimates  is  obtained  from  the  results  of  paper 

number  24when  £^_j(nbn)  ^  <  <®  for  some  q  >  0.  Hence,  the  bandwidth  sequence 

—6 

bR  could  be  of  order  n  for  any  0  <  6  <  1.  This  significant  improvement  for 
general  classes  of  kernel  functions  and  density  functions  (including  multi¬ 
variate  considerations)  is  mainly  due  to  the  probability  cancellation  in  type 
p  spaces  and  the  norm. 

Let  {w  }  be  a  sequence  of  positive  constants  and  W  =  w,+...+w  where 
n  n  1  n 

W  -»  oo  and  w  /  w  -*<*>.  Let  {X  }  be  a  sequence  of  independent  random  elements 
n  n  n  n 

in  D[0,1  J,  In  paper  number  11  the  almost  sure  convergence  of  wjc3Sc 

is  established  under  integral  conditions  and  growth  conditions  on  the  weights 
fwn).  These  results  are  shown  to  be  substantially  stronger  than  the  weighted 
sums  convergence  results  of  Taylor  and  Daffer  (1980),  J ■  Mult .  Anal . ,  and  the 
strong  laws  of  large  numbers  of  Ranga  Rao  (1963),  Theory  of  Probab.  &  Appl.’s, 
and  Daffer  and  Taylor  (1979),  Annals  of  Probab. 
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In  paper  number  26  various  integral  conditions  are  examined  with 
respect  to  the  convergence  almost  surely  and  in  probability  of  the  general 

-CO 

linear  form  /  ^  ank  where  are  random  elements  in  D[0,1  land  {a  is 

an  array  of  real  numbers.  Since  the  hypotheses  of  the  results  of  this  paper 
are  in  general  much  less  restrictive  than  tightness  (nr  convex  tightness), 
these  results  represent  significant  improvements  over  existing  laws  of  large 
numbers  and  convergence  results  for  weighted  sums  of  random  elements  in  D[0,1]. 
Also,  a  counterexample  is  provided  for  both  the  strong  and  weak  laws  of  large 
numbers  when  the  hypotheses  are  weakened. 

Weak  laws  of  large  numbers  in  Banach  spaces  have  typically  been  for 

orthogonal  or  independent  random  elements  and  have  not  been  related  to  the 

geometric  properties  of  the  space.  In  Howell  and  Taylor  (paper  number  6),  a 

class  of  dependent  random  elements  were  studied,  and  weak  laws  of  large 

numbers  for  this  class  (uncondit ional  semibasic  random  elements)  were  obtained 

where  the  convergence  rates  depend  on  the  geometric  structure  of  the  Banach 

rn 

space.  In  particular,  convergence  in  probability  of  ank  Xnfc  *s  obtained 

for  random  elements  {X^}  satisfying  various  distributional  conditions,  including 
independence,  conditional  independence,  and  unconditional  semi-basic,  and 
weights  (a^)  such  that  |a  ^  |  ^  ^  1  for  each  n  where  1  <  p  <  2  is  related 

to  the  geometric  Rademacher  type  of  the  space.  Where  previous  results  required 
type  p  +  6  for  the  space,  these  results  indicate  type  p  is  sufficient.  In 
addition,  these  results  include  a  larger  class  of  random  elements  than  just 
independent  and  conditionally  independent.  Since  lank^  "  *  with  1  <  p  <  2 

leads  naturally  to  applications  to  kernel  density  estimates,  let 

HiT  Ci  "hr5*  ’  Ci  a„k  x„k  ' 

n  n 
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In  particular,  the  use  of  dependent  arrays  of  random  elements  and 

the  convergence  of  weighted  sums  is  more  appropriate  for  the  sequential 

kernel  density est imates  where  the  bandwidth  sequence  may  depend  on  the 

observations  X,,...,X  , ,  that  is,  h  =  h  (X,,...,X 

l  n-  i  n  n  i  n- J 

In  paper  number  7  the  weights  are  again  restricted  to  be  of  the 

1  v  n 

form  *  w^/Wn  and  the  almost  sure  convergence  of  —  £  ^  WaS  °^~ 

n 

tained  under  varying  moment  and  distribution  conditions  on  {x^}.  In  partic¬ 
ular,  if  the  Banach  space  is  of  type  p  and  {X^}  have  uniformly  bounded  tail 
probabilities  by  a  real  random  variable  X,  then 

r  tp_1p[x  >t]  r  dydt<~  (*> 

0  t  yp+l 

is  sufficient  for  the  convergence  where  N(y)  =  card  f  n:  W  /v  <  y\ .  Examination 

n  n 

of  (*)  in  terms  of  moment  conditions  on  {X^}  and  growth  conditions  on  the 

weights  {w^}  was  included.  In  the  absence  of  geometric  conditions  on  the 

Banach  space,  tightness  and  the  integral  condition  in  (*)  are  used  to  obtain 

the  almost  sure  convergence  of  the  linear  forms  when  V  -  w,+. . .+w  .  These 

n  1  n 

results  extend  both  the  real-valued  and  Banach  space  results  of  previous  papers 
by  several  authors. 

C.  Estimation  After  Testing 

In  the  quality  control  of  industrial  or  military  production,  it  is 
common  practice  to  test  hypotheses  to  make  a  decision  about  acceptance  or 
rejection  of  a  production  batch.  Generally,  a  sequential  testing  procedure 
requires  less  data  to  reach  a  decision  than  a  fixed  sample  or  fixed  time  testing 
procedure.  Also,  in  many  situations,  the  experimenter  may  want  to  obtain  a 
point  or  an  interval  estimate  of  some  parameter  after  a  test  of  hypothesis 
has  been  established. 
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Let  Fn  be  a  distribution  function  which  is  stochastically  nondecreasing 

in  0.  In  particular,  F,  can  be  an  exponential  distribution  with  mean  B. 

Several  sequential  procedures  were  developed  which  allow  quick  rejection  of  the 

null  hypothesis  H  :  0  5  9  when  the  alternative  H,  :0s  0,(6,  <  0  )  is  true 
' r  o  o  1  1  1  o 

under  the  general  model  Fg.  In  addition,  these  results  provide  accurate 
interval  estimators  of  the  parameter  0  after  the  test  has  been  completed 
(papers  14,  21). 

D .  Robustness  and  Distribution-free  Procedures 

First,  a  simple  and  robust  estimator  of  the  difference  of  location 
parameters  of  correlated  variables  is  proposed  when  some  observations  on  either 
of  the  variables  are  missing.  We  show  that  this  estimator  is  consistent,  asymp¬ 
totically  normally  distributed,  and  insensitive  to  outlying  observations. 
Asymptotic  relative  efficiency  comparisons  with  other  known  estimators  are 
made  to  show  the  advantage  of  the  proposed  estimator  (paper  9). 

Secondly,  the  robustness  (asymptotic  conservativeness)  of  Kruskal-Wallis 
test  under  certain  departures  from  mutual  independence  of  K  univariate  samples 
is  established  (paper  10).  This  robustness  provides  a  procedure  for  testing 
the  equality  of  K  marginal  distribution  frunctions  based  on  a  broken  random 
sample  from  a  K-variate  distribution  which  satisfied  a  mild  condition.  For 
the  unbroken  sample,  a  generalized  Kruskal-Wallis  test  is  proposed  for  testing 
the  symmetry  of  a  K-dimensional  distribution  function.  The  relative  efficiency 
of  the  K-W  test  against  the  aligned  rank  order  test  is  also  examined  under 
the  normal  shift  model. 

Thirdly,  for  a  randomized  complete  block  design  with  additive  block 
effects,  an  asymptotically  distribution-free  simul t aneous  confidence  region 
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of  pairwise  treatment  differences  is  presented.  The  corresponding  confidence 
bound  has  an  explicit  form  and  is  easily  obtained.  Asymptotic  relative 
efficiency  comparisons  with  other  known  confidence  regions  were  made  to 
show  the  advantage  of  the  proposed  procedure.  An  example  is  provided  for 
illustration  purpose.  The  case  of  treatment  against  control  is  also  discussed 
(papers  27  and  28) . 

Also,  the  problem  of  nonparametric  two-sample  repeated  significance 
tests  was  studied  (paper  number  30).  The  asymptotic  distribution  theory  of 
sequentially  computed  modified  Wilcoxon  scores  was  developed  for  two-sample 
survival  data  with  random  staggered  entry  and  random  loss  tc  follow-up.  A 
repeated  significance  testing  procedure  was  presented  for  testing  the  equality 
of  two  survival  distributions  based  on  the  asymptotic  theory.  The  early 
stopping  properties  of  this  procedure  was  illustrated  by  an  example. 

An  asymptotically  distribution-free  confidence  interval  for  location 
difference  of  two  correlated  variables  with  incomplete  data  on  both  responses 
was  obtained  in  paper  number  17  of  Section  4.  The  confidence  limit  has  an 
explicit  form  and  can  be  easily  obtained.  The  new  procedure  was  also  shown 
to  be  more  robust  than  other  known  parametric  interval  estimators. 

In  paper  number  15  the  point  estimation  of  the  ratio  of  scale  parameters 
in  the  two-sample  problem  with  arbitrarily  right-censored  observations  was  con¬ 
sidered  from  the  minimum  distance  approach.  The  two-sample  Cramer-von  Mises 
distance  was  used.  The  estimator  was  shown  to  be  strongly  consistent  and 
small  sample  properties  wet e  exhibited  by  computer  simulation  results.  Also, 
in  paper  31  nonparametric  point  and  interval  estimators  of  the  ratio  of  two 
scale  parameters  were  given  for  arbitrarily  right-censored  data  based  on  the 
idea  of  Hodges  and  Lehmann.  These  estimators  were  defined  in  terms  of  rank 


15 


test  statistics  for  testing  the  equality  of  two  survival  distributions.  The 
asymptotic  properties  and  efficiencies  of  estimators  corresponding  to  the 
Gehan-Gilbert  test  and  to  the  logrank  test  were  investigated. 

Finally  in  this  area,  paper  number  16  was  a  review  article  on  Friedman’s 
urn  model  written  for  the  Encyclopedia  of  Statistical  Science.  This  urn  model 
has  been  shown  to  be  very  useful  in  the  design  of  sequential  experiments. 

E.  Stochastic  Systems  Theory 

In  paper  number  8  of  Section  4,  general  stochastic  integral  systems 
involving  McShane  stochastic  integrals  were  studied.  These  general  systems 
included  classes  of  stochastic  integro-dif ferential  systems.  Existence, 
uniqueness,  and  stability  of  solutions  to  such  systems  were  investigated 
using  contractor  theory,  which  is  a  generalization  of  the  contraction  principle 
for  nonlinear  operators.  In  particular,  the  results  yielded  frequency-type 
conditions  for  the  stability  of  a  class  of  stochastic  systems. 
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4.  Cumulative  List  of  Written  Publications  in  Technical  Journals 


In  this  section,  the  research  papers  which  have  been  written  under  this 
contract  are  listed.  They  are  divided  into  three  categories:  Published, 
accepted  for  publication,  and  submitted.  Copies  of  the  papers  have  already 
been  forwarded  to  the  Program  Manager  as  they  were  submitted  for  publication 
or  appeared  in  print,  and  sufficient  copies  of  reprints  will  be  forwarded  to 
him  as  the  papers  appear  in  the  future. 


A.  Papers  Published 

1.  W.  J.  Padgett  and  L.  J.  Wei,  Maximum  likelihood  estimation  of  a  dis¬ 
tribution  function  with  increasing  failure  rate  based  on  censored 
observations.  Biomet rika  67  (1980),  470-474. 

2.  K.  F.  Cheng  and  R.  L.  Taylor,  On  the  uniform  complete  convergence  of 
estimates  for  multivariate  density  functions  and  regression  curves. 
Annals,  Inst.  Statist.  Math.,  32  (1980),  187-199. 

3.  L.  J.  Wei,  A  generalized  Gehan  and  Gilbert  test  for  paired  observations 
which  are  subject  to  arbitrary  right  censorship,  Journal  of  the  American 
Statistical  Association  (1980),  634-637. 

4.  L.  J.  Wei,  Recent  development  in  the  design  of  sequential  medical 
trials,  Proc.  Conf.  on  the  Recent  Developments  in  Statistics,  Taipei, 
Taiwan,  (19801,349-356. 

5.  W.  J.  Padgett  and  L.  J.  Wei,  A  Bayesian  nonparametric  estimator  of 
survival  probability  assuming  increasing  failure  rate.  Communications 
in  Statistics,  Theory  &  Methods  A10  (1981),  46-63. 

6.  J.  0.  Howell  and  R.  L.  Taylor,  Marcinkiewiez-Zygmund  type  weak  laws  of 
large  numbers  for  unconditional  random  elements  in  Banach  spaces,  Proc. 
Probab.  in  Banach  Spaces,  111,  LECTURE  NOTES  IN  MATHEMATICS,  Vol.  860 
(1981),  Springer-Verlag,  219-231. 

J.  0.  Howell,  R.  L.  Taylor  and  W.  A.  Woyczynski,  Stability  of  linear 
forms  in  independent  random  variables  in  Banach  spaces,  Proc.  Probab. 
in  Banach  Spaces,  III,  LECTURE  NOTES  IN  MATHEMATICS,  Vol.  860  (1981), 
Springer-Verlag,  231-245. 


7. 
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Papers  Accepted  for  Publication 


A.N.V.  Rao  and  W.  J.  Padgett,  On  the  solution  of  a  class  of  stochastic 
integral  systems,  Journal  of  Integral  Equations  (to  appear  in  1981). 

L.  J.  Wei,  Estimation  of  location  difference  for  fragmentary  samples, 
Biometrika  (to  appear  in  1981). 

L.  J.  Wei,  Asymptotic  conservativeness  and  efficiency  of  Kruskal-Wallis 
test  for  k  dependent  samples,  Journal  of  the  American  Statistical 
Association  (to  appear  in  198 lT^ 

R.  L.  Taylor  and  C.  A.  Calhoun,  On  the  almost  sure  convergence  of 
weighted  sums  of  random  elements  in  DC0,l),  Int.  J.  of  Math  &  Math. 

Scl.  (to  appear,  1981). 

W.  J.  Padgett,  Bayes  estimators  of  reliability  for  the  inverse  Gaussian 
model,  Statistics  T.  R.  /*53 ,  Univ.  of  S.  C.  ,  1980,  IEEE  Transactions 
on  Reliability  (to  appear  in  1981). 

W.  J.  Padgett,  Approximate  prediction  intervals  for  the  mean  of  future 
observations  from  the  inverse  Gaussian  distribution.  Statistics  T.  R. 
#65,  Univ.  of  S.  C. ,  1980,  Journal  of  Statistical  Computation  and 
Simulation  (to  appear) 

W.  J.  Padgett  and  L.  J.  Wei,  A  sequential  test  and  interval  estimation 
in  time  truncated  life  testing,  Statistics  T.  R.  #49,  September,  1979, 
revised  for  Sankhya .  tentatively  accepted. 

W.  J.  Padgett  and  L.  J.  Wei,  Estimation  of  the  ratio  of  scale  parameters 
in  the  two  sample  problem  with  arbitrary  right  censorship,  Statistics 
T.  R.  #69-R,  Univ.  of  S.  C. ,  1980,  under  revision  for  Biometrika. 

L.  J.  Wei,  Friedman's  urn  model,  Encyclopedia  of  Statistical  Sciences 
(to  appear) . 

L.  J.  Wei,  Interval  estimation  on  location  difference  with  inconplete 
data,  Biometrika  (to  appear). 


Papers  Submitted  for  Publication 


K.  Y.  Liang  and  W.  J.  Padgett,  Nonparametric  empirical  Bayes  estimation 
of  reliability.  Statistics  T.  R.  #58,  Univ.  of  S.  C. ,  1980,  submitted  to 
Metron. 

W.  J.  Padgett,  On  estimation  of  reliability  for  the  Bi mbaum-Saunders 
fatigue  life  model,  Statistics  T.  R.  #64,  Univ.  of  S.  C. ,  1980, 
submitted  to  Journal  of  Statistical  Computation  and  Simulation. 
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20.  W.  J.  Padgett  and  L.  J.  Wei,  Conditionally  distribution-free  test  for 
censored  bivariate  observations.  Statistics  T.  R.  #41,  Univ.  of  S.  C. , 
July,  1979,  (submitted  for  publication). 

21.  W.  J.  Padgett  and  L.  J.  Wei,  Interval  estimation  after  sequential 
testing  based  on  the  total  time  on  test.  Statistics  T.  R.  #50,  Univ. 
of  S.  C. ,  October  1979,  submitted  to  Annals  of  the  Institute  of 
Statistical  Mathematics. 

22.  W.  J.  Padgett  and  D.  T.  McNichols,  Estimation  under  reliability  growth 
assuming  gamma  failure  models.  Statistics  T.  R.  #70,  Univ.  of  S.  C. , 
1981,  (submitted  for  publication). 

23.  R.  L.  Taylor,  Convergence  of  weighted  sums  of  random  elements  in  type 
p  spaces,  Statistics  T.  R.  #43,  Univ.  of  S.  C.,  1979,  (submitted  for 
publication) . 

24.  R.  L.  Taylor,  Complete  convergence  for  weighted  sums  of  arrays  of 
random  elements.  Statistics  T.  R.  #48,  Univ.  of  S.  C. ,  1979,  (submitted 
for  publication. 

25.  R.  L.  Taylor,  Tightness  of  measures  and  laws  of  large  numbers  in 
Banach  spaces.  Statistics  T.  R.  #35R,  Univ.  of  S.  C.,  1979,  (submitted 
for  publication). 

26.  R.  L.  Taylor  and  P.  Z.  Daffer,  Stochastic  convergence  of  linear  forms 
in  D[ 0 , 1 ] ,  Statistics  T.  R.  #72,  Univ  of  S.  C. ,  1981,  (submitted  for 
publication) . 

27.  L.  J.  Wei,  Asymptotically  distribution-free  simultaneous  confidence 
region  of  treatment  differences  In  a  randomized  complete  block  design, 
revised  for  the  Journal  of  Royal  Statistical  Society  (B) . 

28.  L.  J.  Wei,  The  analysis  of  a  randomized  complete  block  design  whan 
observations  are  subject  to  arbitrary  right  censorship,  (submitted 
for  publication). 

29.  L.  J.  Wei,  A  conditionally  distribution  free  rank  test  for  testing 
equality  of  distribution  functions  of  correlated  variables  from 
fragmentary  samples,  (submitted  for  publication). 

30.  E.  V.  Slud  and  L.  J.  Wei,  Two-sample  repeated  significance  tests  based 
on  the  modified  Wilcoxon  statistic,  (submitted  for  publication). 

31.  L.  J.  Wei  and  M.  H.  Gail,  Nonparametric  estimation  for  a  scale-change 
model  with  censored  observations,  (submitted  for  publication). 


5. 


Professional  Personnel  Associated  with  the  Research  Effort 


In  addition  to  the  principal  investigators,  Professors  Padgett,  Taylor 
and  Wei,  D.  T.  McNichols  was  a  research  assistant  under  the  contract  and  has 
per formed  the  computing  for  many  examples  used  in  the  research  papers, 
performed  reference  searches,  and  proofread  manuscripts.  This  student  success¬ 
fully  completed  the  written  and  oral  comprehensive  examinations  for  the  Ph.D. 
program  in  September,  1979,  and  began  to  embark  into  an  area  of  research  for 
the  Ph.D.  dissertation  in  the  fall  of  1980,  expecting  to  receive  the  doctorate 
in  the  spring  of  1982. 

Additionally,  during  the  grant  period.  Professor  Taylor  has  directed 
three  students  for  their  degrees  of  Master  of  Science:  (i)  Connie  F.  Sawyer, 
June,  1979,  "The  Markov  Source  as  a  Model  for  Natural  Languages";  (ii)  Patrick 
E.  Flanagan,  May,  1980,  "Computer  programs  for  kernel  estimates  of  a  probability 
density  function";  (iii)  Carol  A.  Calhoun,  May,  1980,  "Laws  cf  large  numbers 
for  D[0,1]  and  estimation  of  density  functions."  Professor  Padgett  has  also 
directed  three  students'  M.S.  theses:  (i)  Y.  S.  Lin  "Confidence  bounds  on 
reliability  for  the  inverse  Gaussian  model,"  completed  in  May,  1980;  (ii)  K. 
Wang,  "Bayes  estimation  of  reliability  for  the  inverse  Gaussian  lifetime  model", 
July,  1980;  and  (iii)  C.  Hinayon,  "Estimation  for  the  Birnbaum-Saunders 
fatigue  life  distribution",  July,  1980.  They  also  were  directing  two  Ph.D. 
students  and  four  additional  masters’  students  in  1980-81,  and  these  students 
will  complete  their  theses  during  the  summer,  1981. 
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6.  Interactions 


The  principal  investigators  presented  (Invited  and  contributed)  papers 
at  various  conferences. 

(i)  W.  J.  Padgett  and  L.  J.  Wei,  "The  maximum  likelihood  estimation  of 
a  distribution  function  with  monotone  failure  rate  based  on  censored 
observations."  The  Annual  Joint  Meetings  of  ASA  and  IMS  in  Washington, 
D.C. 

(ii)  R.  L.  Taylor  (joint  with  P.  Z.  Daffer) ,  "Convergence  of  weighted  sums 

of  random  elements  in  D[0,l3."  Annual  Meeting  of  IMS  in  Washington,  D.C. 

(iii)  L.  J.  Wei,  "A  generalized  Gehan  and  Gilbert  test  for  paired  observations 
which  are  subject  to  arbitrary  right  censorship."  Annual  Meeting  of 
ASA  in  Washington,  D.C. 

(iv)  L.  J.  Wei,  "The  recent  development  of  designs  of  sequential  medical 
trials"  (invited  talk).  Conference  on  Recent  Development  of 
Statistics  at  Academia  Sinica,  Taiwan. 

(v)  W.  J.  Padgett  (joint  with  L.  J.  Wei),  "A  nonparametric  Bayesian  esti¬ 
mator  of  survival  probability  assuming  increasing  failure  rate". 

Eastern  Regional  Meeting  of  IMS,  Charleston,  S.C. 

(vi)  L.  J.  Wei,  "Analysis  of  randomized  complete  block  design  with  arbi¬ 
trary  right  censorship"  (invited  talk).  Eastern  Regional  Meeting  of 
IMS,  Charleston,  S.C. 

(vii)  W.  J,  Padgett  and  Y.  S.  Lin,  "Some  confidence  bounds  on  reliability 
for  the  inverse  Gaussian  model".  The  2nd  Annual  meeting  of  S.C. 

Chapter  of  ASA. 

(viii)  R.  L.  Taylor,  "Complete  convergence  for  weighted  sums  of  arrays  of 
random  elements."  Eastern  Regional  Meeting  of  the  Institute  of 
Mathematical  Statistics,  March  12-14,  1980,  Charleston,  S.C. 

(ix)  R.  L.  Taylor  and  Carol  Calhoun,  "Strong  laws  of  large  numbers  for 
D[0,l3".  The  2nd  Annual  Meeting  of  the  S.C.  Chapter  of  ASA. 

(x)  L.  J.  Wei,  "Asymptotically  distribution-free  simultaneous  confidence 
region  of  treatment  differences  in  a  randomized  complete  block 
design".  The  2nd  Annual  Meeting  of  the  S.C.  Chapter  of  ASA. 

(xi)  W.  J.  Padgett  (joint  with  L.  J.  Wei),  "A  sequential  test  and  interval 
estimation  in  time  truncated  life  testing",  Joint  Statistical 
Meetings,  Houston,  Texas,  August  10-14,  1980. 
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(xii)  R.  L.  Taylor,  "Marcinkiewicz-Zygmund  type  weak  laws  of  large  lumbers 
for  unconditional  random  elements  in  Banach  spaces,"  Invited  Presen¬ 
tation  at  the  Third  International  Conference  on  Probability  in  Banach 
Space,  Tufts  University,  Mass. ,  Aug.  8,  1980. 

(xiii)  W.  J.  Padgett  and  D.  T.  McNichols,  "Estimation  from  development  test 
data  assuming  gamma  lifetimes,"  The  3rd  Annual  Meeting  of  the  S.C. 
Chapter  of  ASA,  April  10,  1981. 

(xlv)  R.  L.  Taylor  and  Carol  Calhoun,  "Weak  laws  of  large  numbers  for 

randomly  weighted  sums  of  random  variables,"  The  3rd  Annual  Meeting 
of  the  S.C.  Chapter  of  ASA,  April  10,  1981. 

(sv)  L.  J.  Wei,  Invited  talk,  Biometrics  Spring  Meeting,  Richmond,  Virginia, 
March,  1981. 

(xvi)  L.  J.  Wei,  Invited  presentation  to  the  Washington,  D.C.,  Chapter  of 
ASA,  May,  1981. 


In  addition,  Padgett  gave  colloquium  talks  at  Auburn  University  and 
the  University  of  Florida,  and  Taylor  gave  a  talk  at  the  University  of  Central 
Florida.  Wei  gave  talks  at  National  Cancer  Institute,  University  of  Maryland- 
College  Park,  University  of  North  Carol ina-Chapel  Hill,  University  of  Wisconsin- 
Madison,  Ohio  State  University,  and  George  Washington  University. 
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7 .  Inventions,  Patent  Disclosures,  and  Applications  Stemming  from  the 

Research  Project 

No  inventions  or  patents  have  stemmed  from  this  research. 

The  results  reported  in  Sections  3  and  4  have  wide  application  in  the 
estimation  and  assessment  of  reliability  of  military  equipment  both  paramet¬ 
rically  and  nonparametrically .  Thus,  maintenance  policies  and  development  of 
new  more  reliable  equipment  may  be  formulated  using  the  statistical  procedures 
and  theory  from  these  general  results.  In  particular,  for  reliability  and  survival 
analysis,  accurate  and  convenient  estimates  of  survival  probabilities,  probability 
density  functions,  and  failure  rate  functions  can  be  obtained  only  assuming  an 
increasing  failure  rate.  Also,  nonparametric  interval  estimation  and  testing 
can  be  performed.  Hence,  assumptions  of  various  parametric  failure  models  are 
not  necessary  for  estimation  or  testing  for  different  types  of  equipment  and 


experimental  situations. 
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8.  Other  Professional  Activities 

During  this  reporting  period,  the  principal  investigators  have  been 
involved  in  numerous  professional  activities  which  are  intimately  related 
to  the  research  efforts  on  this  contract.  The  principal  investigators  have 
refereed  a  total  of  twenty-seven  manuscripts  during  June  1,  1979  to  May  31, 

1981,  three  for  the  Journal  of  Multivariate  Analysis,  one  for  Glasnik  Mate- 

v 

matr icki ,  three  for  IEEE  Transactions  on  Reliability,  five  for  the  Journal 
of  the  American  Statistical  Association,  three  for  Biometrics ,  one  for  the 
Journal  of  Statistical  Inference  and  Planning,  three  for  Communications  in 
Statistics,  three  for  Journal  of  Controlled  Clinical  Trials,  one  for  the 
Canadian  Journal  of  Statistics,  one  for  Technometrics ,  one  for  the  Journal 
of  Econometrics,  one  for  the  Journal  of  Integral  Equations,  and  one  for  the 
International  Journal  of  Mathematics  and  Mathematical  Sciences.  In  addition, 
twelve  reviews  of  papers  were  written  for  Mathematical  Reviews  and  eight  for 
the  Zentralblatt  fur  Mathematik.  Also,  three  NSF  proposals  were  refereed. 

W.  J.  Padgett  and  R.  L.  Taylor  were  the  program  co-chairman  for  the 
Eastern  Regional  Meeting  of  the  Institute  of  Mathematical  Statistics  held  on 
March  12-14,  1980.  They  were  responsible  for  arranging  six  invited  speaker 
sessions  of  the  American  Statistical  Association.  R.  L.  Taylor  was  chairman 
of  a  session  on  "Reliability  and  Growth  Models"  at  the  1980  SREB  Summer 
Research  Conference  on  Statistics  in  Pensacola,  Florida,  and  is  program 
co-chairman  for  the  1982  SREB  Summer  Research  Conference  on  Statistics.  Also, 
Taylor  was  invited  to  talk  at  the  Third  International  Conference  on  Probability 
in  Banach  Spaces  held  August  4-15,  1980,  in  Boston,  Massachusetts.  W.  J.  Padgett 
has  served  as  an  Associate  Editor  of  the  Journal  of  Statistical  Computation  and 
Simulation  since  August,  1980.  L.  J.  Wei  has  served  on  the  Public  Advisory 
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Consnittee  of  NIH  and  chaired  a  session  at  the  1981  Spring  Biometric  Society 
Meeting  in  Richmond,  Virginia. 


